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We investigate the significance of thermai dilepton radiation in the intermediate-mass region 
in heavy-ion reactions at CERN-SpS energies. Within a thermal fireball model for the space-time 
evolution, the radiation from hot matter is found to dominate over hard 'background' processes 
(Drell-Yan and open charm) up to invariant masses of about 2 GeV, with a rather moderate fraction 
emerging from early stages with temperatures T ~ 175—200 MeV associated with deconfined matter. 
Further including a schematic acceptance for the NA50 experiment we find good agreement with 
the observed enhancement in the region 1.5 GeV < M MM < 3 GeV. In particular, there is no need 
to invoke any anomalous open charm enhancement. 
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The identification of experimental signals from a possibly formed Quark-Gluon Plasma (QGP) constitutes one of 
04 , the prime goals in the (ultra-) relativistic heavy-ion program. Among the earliest suggestions [[[] for QGP signatures 
has been the appearance of an enhanced radiation of dileptons (and/or photons) in the intermediate-mass region 
(IMR), 1.5 GeV < Mu < 3 GeV, i.e., well above the light vector mesons p, lu and <f> but below the J/\& resonances. 
The importance of this region resides on the facts that, on the one hand, low-energy hadronic processes (such as 
7T7T annihilation or Dalitz decays) are sufficiently suppressed, and, on the other hand, hard processes (in particular 
• Drell-Yan annihilation), which prevail in the high- mass region M > 4 GeV, increase rather slowly towards smaller M 
Q\ ' and thus may be over-shined by thermal radiation. 

At the CERN-SpS, intermediate-mass dilepton spectra have been measured in the dimuon mode by the 
G\ ' NA38/NA50 |||| and HELIOS-3 || collaborations (the dielectron data from the CERES experiment § lack sufficient 
statistics in this region). In central collisions of heavy nuclei both experiments found factors of 2-3 enhancement over 
■ the extrapolation of known sources from proton-induced collisions, given by primordial Drell-Yan annihilation || as 
Oh! wen as semileptonic decays of associatedly produced charmed (D, D) mesons. 

Several suggestions to explain this excess have been made, e.g., an enhanced production of cc pairs. However, such an 
| increase is not easily justified theoretically , and would also be difficult to reconcile with the current understanding 
of the observed J/^> suppression. Another possibility has been pursued by Lin and Wang || who showed that 
Z?-meson rescattering in hot/dense matter (without invoking anomalous enhancement) might generate a transverse 
• i-h , momentum broadening which can enrich the p + p~ phase space covered by the NA50 experiment. The resulting 
increase is, however, much too small to explain the data ||. Spieles et al. ^ have investigated the role of secondary 
Drell-Yan processes in hadronic rescatterings within the UrQMD transport model and found that at an invariant 
mass of, e.g., = 1.5 GeV this contribution may constitute up to ~ 30 % of the primordial Drell-Yan yield. This 
is again far from being able to account for the experimental findings. Finally, Li and Gale analyzed the IMR of the 
HELIOS-3 data using a transport model incorporating dilepton production through secondary hadronic annihilation 
processes p0[ . They found that the enhancement in S+W collisions can indeed be explained; unfortunately, the 
HELIOS-3 data quickly run out of statistics for > 2 GeV. Since the transport framework employed was entirely 
based on hadronic degrees of freedom, a possible Quark-Gluon Plasma formation was not addressed explicitly. 

Another important issue is the consistency with dilepton measurements in the low-mass region (LMR). The 
CERES/NA45 S and HELIOS-3 Q collaborations have reported a significant excess over final-state hadron de- 
cays which can indeed be explained by thermal radiation once medium effects in the vector meson spectral functions 
(especially for the p) are accounted for |Ilj , |l2"| (see Ref. (Rj for a recent review). In this letter we would like to 
investigate whether thermal production in the 1.5-3 GeV region, folded over the same fireball evolution as employed 
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for the LMR in Ref. [ 
formation, we study t 



2[ , can simultaneously describe the NA50 data. In addition, by allowing for a possible QGP 
re significance of dileptons produced through quark-antiquark annihilation in the early phases 



of central Pb(158 AGeV)+Pb collisions. 
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II. DILEPTON SOURCES AT INTERMEDIATE MASSES 



In the absence of any medium (or rescattering) effects, dilepton production in hadronic collisions arises from either 
initial hard processes or the decay of produced mesons. At invariant masses beyond 1.5 GeV, the contribution from 
light hadron decays (containing u, d and s quarks) sharply drops off either due their insufficient mass (as, e.g., for 
p, to and <f> mesons) or because of their increasing total widths resulting in a small electromagnetic branching ratio 
(as, e.g., for the p 1 resonances). For masses between 1.5 GeV < M < 3 GeV the most important decay contribution 
therefore arises from the associated decay of charmed mesons D, D originating from primordial production of a cc 
pair (the thermal abundance of D-mesons is negligible at SpS energies). The other significant 'background' source in 
the IMR is the well-known Drell-Yan (DY) process H . The latter entirely saturates the yield at high masses bey ond 
~5 GeV and can thus be exploited for normalization purposes as has been verified by the NA50 collaboration p4|| . To 
be able to compare our calculations to the NA50 data we will adopt the same strategy here (in addition, the shape of 
the DY-yield after applying acceptance cuts can be checked versus experimental simulations using event generators). 
In a central collision of two equal nuclei with mass number A the number of produced Drell-Yan pairs per invariant 
mass and rapidity interval can be obtained as 

DY (h = 0\= ,-1 4/3 DY m 
dMdy K ' \-kR\ dMdy 1 ' 

in terms of the (standard) elementary Drell-Yan cross section in a nucleon-nucleon collision, 
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Here, q{x\^ 2 ) (9(^1.2)) denote the (anti-) quark distribution functions within a pair of colliding nucleons from the 
two incoming nuclei (neglecting nuclear effects). Their arguments are related to the center-of-mass rapidity y and 
invariant mass of the lepton pair as xi y2 = xe ±y with x = M/^fs, where s denotes the total cms energy of the 
nucleon-nucleon collision. For the parton distribution functions we employ the recent leading order parameterization 
from Ref. (=GRV-94 LO), which necessitates a A'-factor of 1.5 [pf (accounting for higher order effects) to comply 
with p-A Drell-Yan data. As has been further pointed out in Ref. [p|, the explicit inclusion of isospin asymmetric 
sea-quark distributions in the GRV-94 LO set entails substantially smaller isospin corrections for nuclear collisions; 
e.g., for Pb+Pb the correction is as small as 3 %, which will be neglected in the following. In Eq. (|l]), the root-mean 
squared radius parameter Rq ~ 1.05 fm arises from a folding over a Gaussian thickness function Pq| . For (slightly) 
non-central collisions as considered below we simply replace A in Eq. ([[J) by the number of participants. 

In nuclear collisions additional so-called 'thermal' dileptons are radiated through interactions during the expan- 
sion/cooling of the hadronic fireball. In thermal equilibrium, the dilepton production rate depends on the temperature 
and (baryon-) density of the matter and is governed by the electromagnetic current-current correlator II em . Whereas 
at low masses medium modifications of the light vector mesons (especially the p) lead to a strong reshaping of the 
correlator [ ^2| , the situation simplifies in the intermediate-mass region: already in vacuum the cross section for the 
reverse process of e + e _ — > hadrons is well-described by lowest-order perturbative annihilation into qq pairs; at the 
correspondingly probed distance scales of d < l/(1.5GeV) ~ 0.13 fm, additional medium effects are expected to be no 
longer relevant. In fact, as has been shown in Ref. []l0|| , a rate calculation involving a large set of hadronic reactions 
fitted to reproduce exclusive e + e~ — > hadrons cross sections yields an equivalent dilepton production rate (within the 
experimental uncertainties), i.e., the hadronic and partonic description of II em are 'dual' for M > 1.5 GeV. We will 
thus employ the simple partonic rate 
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q—u,d,s 

= ^f B (qo;T) 2 - (3) 

and note that in the IMR both temperature and density corrections, being of order 0{T/M) and 0(p q /M), are 
irrelevant under the conditions probed at the SpS. To obtain the total yield from thermal radiation, the rate (||) 
has to be integrated over the space-time history of a given collision system. Here we employ a simplified thermal 
fireball description including, however, the most important features of a full hydrodynamic simulation in a transparent 
way. As delineated in Ref. |l2] ] we first construct a trajectory in the baryon-chemical potential/temperature plane 
based on baryon-number and entropy conservation. Using a resonance gas equation of state (EoS) including the 32 
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(16) lowest lying mesonic (baryonic) states with s/gs=26 it is found to be consistent with the recently determined 
chemical freezeout point of (hb, T) c hem = (270 ± 20, 175 ± 10) MeV (Tt]] , which we also assume to be the critical point 
separating the hadronic and quark-gluon phase. In the hadronic evolution towards thermal freezeout, (T, fix) f = 
(115 ± 10, 430 ± 30) MeV, one additionally requires pion (kaon) number conservation to maintain the final pion-to- 
baryon ratio of 5:1 at SpS energies, entailing the build-up of finite pion (kaon) chemical potentials reaching fj,„ = 
60 — 80 MeV (/j, k ^ = 100 — 130 MeV). For the initial stages we assume an ideal QGP EoS. A time scale is introduced 
by employing a cylindrical hydro-type volume expansion 

V FC = 2{z + v z t + \a z t 2 ) 7T (r + ^ ± t 2 ) 2 (4) 

with a forward and a backward going firecylinder |fl2f to accomodate a realistic rapidity spread in the final particle 
distributions. r$ corresponds to the initial transverse overlap of the two colliding nuclei (determined by the impact 
parameter b), and the parameters v z , a z and a± are adjusted to the finally observed flow velocities. Given the total 
entropy of the reaction, the value of Zq is equivalent to a formation time (or initial temperature). Once the system has 
cooled down to T c , the fraction / of matter in the hadronic phase can be inferred from the standard entropy balance 

fs HG {Tc) + (1 - f)s QG p(T c ) = Stot/V FC (t) . (5) 
The time-integrated dilepton yield from one collision event then takes the form 

tfo OO 

^jytherm Ct 2 f fa 2 

—777- = o— M / dt V Fc(t) / dq *- f B (q ;T(t),^(t)) Acc(M,g,y) , (6) 
dM 3 7r J J J q 

o o 

where the factor Acc(M, q, y) simulates the experimental acceptance that the produced muons undergo in the NA50 
detector, see below. In the hadronic phase, the Bose factor f B (qo] T(t), n^{t)) includes the enhancement due to finite 
/Utt's. Since in the IMR the vector correlator is dominated by nai annihilation Jl(| , which is a 4-pion type process, 
the ^ = value of the Bose factor is augmented by a factor ~ exp(4/i7r/T). 

Before we come to the full results let us compare the relative importance of Drell-Yan versus thermal pairs without 
any cuts. Fig. [j] shows the total dilepton yields from Drell-Yan annihilation and the time-integrated thermal signal 
from a b — 1 fm Pb(158 AGeV)+Pb collision. The latter has been calculated starting from an initial temperature 
Ti = 192 MeV (corresponding to zo = 2.3 fm in Eq. (H)) cooling down to about T/ Q = 113 MeV within a total fireball 
lifetime of ~ 14 fm/c. As expected, the thermal curve (solid line) dominates towards low invariant masses but drops 
below the DY component (dashed-dotted line) at M > 1.8 GeV due its much softer slope. 




0.0 1.0 2.0 3.0 4.0 5.0 6.0 
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FIG. 1. Drell-Yan (dashed-dotted line) versus thermally produced fi + /j, pairs in central Pb(158AGeV)+Pb, integrated over 
all rapidities and without any acceptance cuts. 
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To test the sensitivity with respect to the initial conditions, we have carried out a calculation with zq = 1.5 fm which 
translates into an initial temperature Tj = 221 MeV. With the same fireball lifetime of 14 fm/c as before (leading to 
a practically identical freezeout temperature) the thermal spectrum looks rather similar to the Ti = 192 MeV case 
up to M ~ 2 GeV, indicating that the decrease of the rate at smaller temperatures (corresponding to later stages) is 
essentially compensated by the larger volumes. Only at masses beyond 2 GeV the higher initial temperature entails 
a visible increase of the thermal yield together with a slightly harder slope. 



III. IMR SPECTRA IN NA50 



In order to compare our calculations to the measurements of NA50 we have to account for the experimental 
acceptance cuts. In the IMR the geometry of the spectrometer imposes kinematic restrictions which are approximately 
given by |J) 

o< y c ™ s <1 

-0.5 < cosQ+ s < 0.5 , (7) 
where y c ™ is the pair rapidity in the cms of the colliding nuclei and 

cos9+ = P % Pz . (8) 

[(1-/T ir )(MV4-m2)]i/2 

Here, Q^S denotes the polar angle of the /i + with respect to the z-axis in the rest frame of the dimuon (the so-called 
Collins-Soper frame) with pf, E + and — (pf + p~)/(E + + E~) measured in the laboratory frame. Furthermore, 
the single-muon tracks in the lab-frame are selected according to 

0.037 < 9^ < 0.108 . (9) 

More difficult to asses are absorption effects on the single-muon tracks; as a very rough characterization thereof we 
require the single-muon energies to satisfy [ fl8| 

16000(9^ - 0.065) 2 , 0.037 < 9 M < 0.065 

,0.065 < 9 P < 0.090 (10) 
13000(9^ - 0.090) 2 , 0.090 < 9 P < 0.108 , 




where E cut has been chosen as 8 GeV in Ref. [|3|. Following Refs. Jl 9|j3|] , we have also incorporated a phenomenological 
qt-broadening of the Drell-Yan spectra by supplying Eq. ([!]) with a Gaussian factor exp(— / 2a^ t ) / 2a^ t to obtain 
dNpY I 'dM dydq^ with a qt — 0.8 — 0.9 GeV. Finally, we account for the finite experimental mass resolution by folding 
the calculated spectra with a Gaussian, the effective <jm being estimated from the apparent width of the J/^ peak. 



NA50 

sim. (NA50) 
E M =9GeV 
E cu =12GeV 
E„„=15GeV 
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FIG. 2. Effects of the schematic cuts, Eqs. (Q) and (|lo|), on the invariant mass spectrum of Drell-Yan produced pairs in 
central Pb(158AGeV)+Pb. In comparison to a detector simulation of the NA50 collaboration (dotted line; cf. also Figure 1 in 
Ref. or the Drell-Yan fits shown in Ref. ||)) the curve calculated with E cu t — 12 GeV seems to give the best representation 
of the experimental acceptance. 
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We have checked the schematic acceptance by comparing our results for the Drell-Yan production with and without 
the cuts to a simulation of the NA50 collaboration using the PYTHIA event generator M] . The pertinent results are 
displayed in Fig. || using three different values for the single-muon energy threshold. The latter does not significantly 
affect the high-mass tail beyond 4 GeV which has been used to fix the normalization in comparison to the NA50 
data. From the behavior towards lower masses it seems that E cut ~ 12 GeV gives the closest reproduction of the 
simulation results from the NA50 collaboration 19 gj. We will therefore use this value in the following also for the 
thermal radiation contribution, but note that this might not be an accurate representation of the actual experimental 
situation. 

To enable a more complete comparison with the NA50 data, the associated decays of open charm mesons have to 
be included. We simply adopt this contribution as given in Ref. jj| (without any anomalous enhancement factor). 

In Fig. |^ our final results are confronted with the NA50 data in central Pb+Pb collisions at 158 AGeV projectile 
energy. Given the uncertainties in the simulation of the experimental acceptance, the excess observed in the IMR 
(1.5 GeV < M < 2.5 GeV) can be well accounted for by thermal radiation when combined with Drell-Yan and open 
charm pairs (resulting in the full line in the left panel). Also the shape and magnitude of the transverse momentum 
spectra (right panel of Fig. |§|) are rather well described once thermal pairs are included. 




FIG. 3. IMR dilepton spectra in comparison to NA50 data from central Pb(158 AGeV)+Pb collisions with N part — 380 (left 
panel: invariant mass spectra, right panel: transverse momentum spectra) . The short-dashed and dashed-dotted curves are the 
calculated thermal and Drell-Yan yields, respectively, using the approximate acceptance as given in the text. The long-dashed 
curve in the left panel is the part of the thermal contribution originating from the plasma phase. The dotted curves represent the 
open charm contribution as obtained by the NA50 collaboration in a PYTHIA simulation without anomalous enhancement Q . 
The full lines are the sum of thermal, open charm and Drell-Yan dileptons. 



As mentioned in the introduction an important issue in the IMR is the quantitative role of radiation from the 
QGP. With our standard value for the initial temperature of Ti — 192 MeV, the qq contribution to the thermal yield 
constitutes about 25 % (slightly more if a higher Ti is assumed) . This is significantly more than typical estimates for 
the low-mass region, but still moderate owing to the rather small space-time volume occupied by a QGP at the (full) 
SpS energies. Although perturbative a s -corrections might enhance the partonic production rate somewhat the major 
signal will still originate from the hadronic phase. 

One may further ask to what extent the presence of finite pion chemical potentials influences the thermal yield. 
When using fj, n = throughout the latter is reduced by about 30-40 % (we should stress, however, that such 
a calculation underestimates the measured pion multiplicity in central Pb(158 AGeV)+Pb by ~ 40-50 %). Thus, 
similar to what has been found in the LMR Jl2| , chemical off-equilibrium effects (determined by imposing pion number 
conservation) seem to be an important ingredient to understand dilepton production in the IMR. 



IV. CONCLUSIONS 



In summary, we have evaluated dilepton spectra in central Pb(158 AGeV)+Pb collisions as measured by the NA50 
experiment. At intermediate masses the total spectra are composed of 3 major components: open charm decays, 
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Drell-Yan annihilation and thermal radiation. Whereas the first one has been taken from NA50 event generator 
simulations, the Drell-Yan part has been calculated explicitly to both determine the data normalization and check 
the approximate acceptance of the NA50 detector. The thermal contribution is based on a realistic fireball expansion 
in accord with recent hadro-chemical analysis. Our main result is that the excess observed by the NA50 collaboration 
in the mass region 1.5 GeV < M < 2.5 GeV can be explained by the thermal signal without invoking any anomalous 
enhancement in the charm production. The importance of the thermal yield quickly ceases beyond M ~ 2.5 GeV, thus 
not affecting the interpretation of J/^> suppression effects - as opposed to scenarios invoking open charm enhancement. 
The contribution from the QGP phase turned out to be around 20-30 %. 

Since the low-mass dilepton spectra from CERES /NA45 can be explained in the same framework (including medium 
effects in the low-lying vector mesons) a consistent picture of dilepton production at the full CERN-SpS energy seems 
to emerge. Similar conclusions have been reached in Ref. pCfl where a more schematic dynamical evolution has been 
employed. 

Let us finally address the question what can be done from the experimental side to further test our results. Firstly, 
a direct measurement of D mesons will fix the open charm contribution to the dilepton spectra. Secondly, a dedicated 
(short) run of the NA50 experiment with reduced magnetic field to focus on the IMR would be most valuable. The 
resulting gain in acceptance for low-qt dileptons might quantitatively disentangle the thermal signal from hot/dense 
matter, putting model predictions under scrutiny and increasing the sensitivity to the very early stages of the collision. 
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